Anomalous Meissner effect in NS junction with spin-active interface 
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We investigate Meissner effect in normal metal/superconductor junctions where the interface is 
spin-active. We find that orbital magnetic susceptibility of the normal metal shows highly nontrivial 
behaviors. In particular, the magnetic susceptibility depends on the temperature in an oscillatory 
fashion, accompanied by its sign change. Correspondingly, magnetic field and current density can 
spatially oscillate in the normal metal. The possible spontaneous formation of the current pattern 
is also discussed. These results are attributed to the generation of odd- frequency pairing due to the 
spin-active interface. 
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Interface phenomena related to the superconductivity 
constitute a rich field of condensed matter physics. When 
superconductor is attached to normal metal, Cooper 
pairs penetrate into the normal metal which acquires su- 
perconducting correlation. This is called the proximity 
effect. As a result, for example, the normal metal has a 
gap in the density of states or shows Meissner effect P-|5|. 
In most cases, when lowering temperature, the proximity 
effect and hence the Meissner response become stronger. 

However, unexpected behavoir of the proximity in- 
duced Meissner response has been reported: Mota et al. 
have uncovered a low-temperature anomaly in the mag- 
netic response of cylindrical structures. At very low tem- 
peratures, the susceptibility shows a reentrant behavior 
and even has paramagnetic region. [6- 12] However, the 
origin of this phenomenon still remains unclear. 

Recently, it has been clarified that in normal 
metal/superconductor junctions, if the interface is spin- 
active, induced superconducting pairing in the normal re- 
gion can change its symmetry from even-frequency pair- 
ing to odd- frequency pairing. Q Here, even- or odd- 
frequency means that Cooper pair wavefunction is even 
or odd with respect to Matsubara frequency (or imag- 
inary time).[l5| If proximity induced pairing symmetry 
changes, the associated Meissner effect will also change 
qualitatively. This is the problem we address in this pa- 
per. 

In this paper, we study Meissner response in the nor- 
mal metal attached to superconductor where the inter- 
face is spin-active. We find that orbital magnetic suscep- 
tibility of the normal metal shows quite complex depen- 
dence on junction parameters. In particular, the mag- 
netic susceptibility depends on the temperature in an 
oscillatory fashion, accompanied by its sign change. We 
also show the behavior of magnetic field and current den- 
sity. These can spatially oscillate in the normal metal. 
These results are attributed to the generation of odd- 
frequency pairing which stems from the spin-active in- 
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FIG. 1: (Color online) Schematic picture of the model. Left: 
arrow in the circle represents the direction of spin, which is 
rotated at the scattering by the spin-active interface. Right: 
structure of the spin- active interface at x = L. 

We consider a junction consisting of a diffusive nor- 
mal metal (DN) with a length L and resistance R^^ and 
a superconductor. We schematically show the model in 
Fig. [TJ The interface between the DN and the super- 
conductor at X = L has a resistance Rh (or tunneling 
conductance Gt) and the surface at x = is specular. A 
weak external magnetic field H is applied in z-direction 
(perpendicular to the plane of left panel of Fig. 1). We 
consider spin-active interface at x = L which is described 
by mixing conductances and which reflect the spin 
rotation upon reflection and transmission at the inter- 
face, respectively. ^3 To evaluate and G^, we 
model magnetic barrier (interface) region as a rectangu- 
lar potential V with the exchange field h and the width d 
as shown in the right panel of Fig. [TJ following Ref.[l3]- 
Figure [2] shows mixing conductances (upper) G^^IGt and 
(lower) G^/Gt as functions of kpd and V/Ep with the 
Fermi wave vector of the DN kp and the Fermi energy Ep- 
Mixing conductances oscillate with these parameters and 
G(f)/GT rapidly increases when approaching V/Ep — 1. 

To study the Meissner response, we adopt the qua- 
siclassical Green's function theory. The normal and 
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FIG. 2: (Color) Mixing conductances (upper) G^IGt and 
(lower) G^IGt for hjEp = 0.01. 



anomalous Green functions are parameterized as = 
cosher and = singer, respectively, where cr(= ± =t7i) 
denotes spin. The Usadel equation [18] in the DN reads 



D—^9a - 2uJn sin6>a 







(1) 



where D and uo^ are diffusion constant and Matsubara 
frequency, respectively. The boundary conditions are de- 
termined by the continuity of the matrix current, and are 
given by 

_a_ 
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at X = L[17\. Here, g and / are bulk Green's functions in 
the superconductor. Mixing conductances appear when 
spin symmetry is broken in the interface region. [17] Con- 
sequently, spin triplet pairing is generated, which should 
be odd in frequency in DN because only 5-wave pairing 
can survive impurity scattering, which is the basic as- 
sumption of the Usadel equation. [13] 

When a magnetic field is applied parallel to the inter- 
face, rich and nontrivial screening effect occurs. Within 
the linear response theory, the current distribution flow- 
ing in ^/-direction is given by [H, Q 



-S7Te^N{EF)DT ^ sin^ ( 



{x)A{x), (4) 



where A(x), N{Ef) and T denote the vector potential, 
the density of states at the Fermi energy and the temper- 
ature of the system, respectively. The Maxwell equation 
reads 



dx'^ 



A {x) = -Anj {x) . 



The boundary conditions for A{x) are given by 



d 

dx 



A (0) = H, A (L) = 0, 



(5) 



(6) 



where we have neglected the penetration of magnetic 
fields into the superconductor by assuming a small pen- 
etration depth in superconductor. 

Finally, we obtain the expression of the orbital mag- 
netic susceptibility. 



47rx = 1 



A{0) 



(7) 
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We set h/Ep = 0.01, Rd/Rb = 10 and 
IGne^N (Ep) = 1000. Below, ^ and Tc denote the su- 
perconducting coherence length and the transition tem- 
perature, respectively. In the following, we plot — 47rx 
with its magnitude less than unity since |47rx| > 1 state 
indicates an instability toward some ordering or sublin- 
ear dependence on magnetic field due to the breakdown 
of linear response theory. When — 47rx > 1, the per- 
meability /i = 1 + 47rx becomes negative and hence the 
energy density 5^/2/i is unstable at 5 = 0. This sug- 
gests that the new ground state with spontaneous current 
and magnetic field distribution is stabilized. Also, cor- 
respondingly, plots of magnetic field and current density 
are restricted to a certain regime of magnitude. 

Figure [3] shows susceptibility — 47rx at L/^ = 10 with 
(upper) V/Ef = 0.95 and (lower) V/Ef = 1. At 
V/Ef = 0.95, over some region, paramagnetic state, 
namely that with positive x^ appears. At V/Ef = 1, 
stronger oscillation of the susceptibility accompanied by 
its sign change is seen. When mixing conductance is 
present, odd- frequency pairing is generated in the DN 
region. 13[ [l^ This odd-frequency pairing makes it pos- 
sible to oscillate magnetic field rather than suppress in 
the DN region [l9|. as expicitly shown below. Therefore, 
susceptibility could be positive when odd-frequency pair- 
ing correlation is dominant over even-frequency pairing 
in the DN. 

If purely even (odd) frequency pairing state is realized 
in the DN, 6a- becomes purely real (imaginarv).[l9[ [2q| 
Then, the sign of screening current Eq.(j4]) depends on 
whether induced pairing in the DN is even- or odd- 
frequency pairing. This drastically changes the suscep- 
tibility. To understand this qualitatively, let us consider 
thin limit of DN where spatial dependence of Oa is negli- 
gible. Then, for purely even-frequency pairing state, the 
Maxwell equation reads 



dx"^ 



A{x) = k^A{x) 



(8) 
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FIG. 3: (Color) Susceptibility -Attx at L/^ = 10. (upper) 
V/Ef = 0.95 and (lower) V/Ef = 1. 



with a real constant k. Then, we have 



H 



A{0) = — - tdnih kL. 
k 



(9) 



Upon insertion of this equation into Eq.([7|), we find that 
— 47rx is positive definite. On the other hand, for purely 
odd- frequency pairing state, the Maxwell equation reads 



with a real constant Hi. Then, we obtain 

A{0) tan/^L. 

hi 



(10) 



(11) 



Substituting this equation into Eq.Q, we find that — 47rx 
can change its sign. Moreover, it can show divergent be- 
havior near i^L = 7r/2 mod tt. In this way, we can un- 
derstand the behavior of the susceptibility in Figure [3l 
Positive sign of x means that proximity induced super- 
conductivity shows paramagnetism. It has been known 
that in d-wave superconductors, paramagnetic contribu- 
tion to Meissner effect arises from the Andreev surface 
states. |2ll-l24| In stark contrast, paramagnetic Meiss- 



ner effect predicted here does not require unconventional 
superconductivity. 

Figure H] shows susceptibility at kpd = 10 with (up- 
per) V/Ef = 0.95 and (lower) V/Ep = 1. The suscepti- 
bility oscillates with temperature while L dependence is 



FIG. 4: (Color) Susceptibility -47rx at /cfo? = 10. (upper) 
V/Ef — 0.95 and (lower) V/Ef — 1- The inset shows sus- 
ceptibility as a function oi T/Tc for L/^ — kFd = 10 and 
V/Ef = 0.95. 



weaker. The inset shows susceptibility as a function of 
T/Tc for L/^ = kpd = 10 and V/Ep = 0.95. Refiect- 
ing the presence of odd-frequency superconductivity, the 
susceptibility shows oscillating divergent behavior (The 
border between red and green regions in the main panel 
corresponds to the divergence of the susceptibility). 

Figure [5] depicts normalized magnetic field H{x)/H 
(upper) and current density (lower) at L/^ = kpd = 10 
and V/Ef = 0.95. The current density j{x) is plotted 
in the unit of TcLH/2D. For Tc ImeV, L - 1/im, 

- 10 and - lO-^m/s, TcLH/2D - lO^^A/m^. As 
seen, both magnetic field and current oscillate in space. 
In a similar way to obtain Eq.(ll), we can show that the 
oscillation is due to the odd-frequency pairing. There- 
fore, Figure [5] also indicates that odd- frequency pairing 
does not repel magnetic field. The relation between sus- 
ceptibility and magnetic field can be obtained along the 
same line: 



•47rx 



tan tvL 
nL 



1 



sin nL H{L) 



(12) 



This indicates that, to realize paramagnetic state, H{L) 
has to be larger than i7, namely, magnetic field at the 
interface should be larger than that at the surface of the 
DN. 

To attain anomalous Meissner effect, dominant odd- 
frequency pairing in the DN is required. This can be 
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FIG. 5: (Color) Magnetic field (upper) and current density 
(lower) at L/^ = kpd = 10 and V/Ef = 0.95. The plots 
are restricted to a certain regime of magnitude for clarity of 
figure. 

acheived when magnitude of mixing conductance is com- 
parable to that of tunneling conductance (Gt) ll3|]. As 
seen from Figure [21 to enhance effect of spin-activeness, 
low interface transparency, namely, large potential bar- 
rier V is desirable. This may be acheived by using ferro- 
magnetic insulator such as EuO, EuS or La2BaCu05 as 
a tunneling barrier. 

We find divergent behaviors of the susceptibility, the 
magnetic field and the current density. However, their 
magnitudes would be reduced due to sgin orbit scattering 
or magnetic scattering in the sample [25] but the behav- 
iors presented above are robust against this effect within 
the reasonable range of 1/{TcTs) < 1 with being the 
lifetime of the spin. 

Finally, stability and Meissner effect of odd-frequency 
pairing in the bulk have been a controversial issue (see 
Refs. [26|, l27j). However, the presence of odd- frequency 
pairing correlation in this paper is due to spin symme- 
try breaking at the interface and thus irrelevant to these 
problems. Namely, the proximity effect is uniquely de- 
termined by the boundary conditions. 

In summary, we have studied Meissner response in the 
normal metal attached to superconductor where the in- 
terface is spin-active. We found that due to the spin- 
active interface, the magnetic susceptibility shows very 
complex behaviors. In particular, the susceptibility de- 
pends on the temperature in an oscillatory fashion, ac- 
companied by its sign change. We also showed that mag- 
netic field and current density can spatially oscillate in 
the normal metal. The possible spontaneous formation 



of the current pattern was also discussed. These results 
are attributed to the generation of odd-frequency pairing 
arising from the spin-active interface. Our results could 
be confirmed by experiments with /i-SR or microwave 
resonance. 
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